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LONG-TERM  GOALS 

The  long  term  goal  is  to  develop  the  capability  to  predict  the  internal  wave  climate  in  the  ocean  with 
the  high-resolution,  three-dimensional,  nonhydrostratic  SUNTANS  model. 

OBJECTIVES 

Our  primary  objective  is  to  understand  the  mechanisms  leading  to  the  generation  of  internal  waves  in 
the  South  China  Sea  by  performing  high-resolution,  nonhydrostatic  numerical  simulations  with  the 
SUNTANS  model.  A  secondary  objective  is  to  perform  simulations  to  aid  with  the  planning  of  the 
field  experiments  for  the  IWISE  (Internal  Waves  in  Straits  Experiment)  DRI. 

APPROACH 

We  are  employing  a  combination  of  theoretical  analyses  and  numerical  simulations  with  the 
SUNTANS  (Fringer  et  al.,  2006)  model  to  perform  a  detailed,  nonhydrostatic  study  of  the  mechanisms 
governing  the  generation  of  internal  waves  in  the  Luzon  Strait,  which  forms  the  eastern  boundary  of 
the  South  China  Sea  (see  Figure  1).  Our  work  builds  upon  simulations  we  have  performed  as  part  of 
the  NLIWI  (Nonlinear  Wave  Initiative)  DRI  that  have  been  designed  to  understand  the  generation  and 
propagation  of  internal  waves  in  the  South  China  Sea,  as  well  as  how  the  internal  waves  interact  with 
the  continental  shelf  on  the  western  side  of  the  sea.  The  simulations  are  similar  to  the  NLIWI 
simulations  but  focus  on  the  Luzon  Strait.  We  plan  on  performing  simulations  in  smaller-scale 
domains  in  order  to  focus  on  the  generation  region,  rather  than  full  basin-scale  domains,  which  are 
better-suited  to  studying  internal  wave  generation  and  the  deep-basin  propagation  and  interaction  with 
the  Dongsha  Plateau  and  the  Chinese  continental  shelf  (see  Figure  1).  SUNTANS  domains  for  the 
present  work  are  targeted  towards  simulating  high-resolution,  nonhydrostatic  flow  physics  through 
channels  and  over  sills  to  understand  the  fundamental  physics  of  internal  wave  generation.  We  expect 
the  grid  to  be  refined  around  the  Babuyan  and  Batan  Islands  in  order  to  highly  resolve  the  tidal  flow 
around  these  islands  and  to  resolve  possible  hydraulic  control  features.  Based  on  our  experience  with 
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simulating  internal  waves  in  the  SCS,  we  expect  to  be  able  to  use  SUNTANS  in  a  predictive  manner  to 
aid  in  the  placement  of  instrumentation  for  the  field  component  of  the  IWISE  DRI. 
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Figure  1:  Bathymetric  map  of  the  South  China  Sea,  showing  the  locations  of  mooring  sites  LI, 
B1  and  S7  from  the  WISE/VANS  project  in  2005,  and  the  transects  that  are  used  for  the  two- 

dimensional  theoretical  analysis. 


WORK  COMPLETED 

We  have  developed  a  linear  theory  for  the  generation  of  internal  waves  over  ridges  to  determine  the 
phase  in  the  tide  at  which  peaks  in  the  depression  waves  cross  the  ridges. 

RESULTS 

The  most  basic  aspect  of  nonlinear  internal  waves  in  the  South  China  Sea  that  has  yet  to  be  explained 
is  the  correlation  of  peaks  in  westward-propagating  internal  waves  with  maximum  eastward  (ebb) 
currents  at  the  eastern  ridge  in  the  Luzon  Strait.  Ramp  et  al.  (2004)  used  linear  theory  to  compute  the 
internal  wave  phase  speeds  and  trace  peaks  in  the  internal  wave  signatures  observed  at  B 1  back  to 
mooring  LI,  and  showed  that  the  tracebacks  always  crossed  mooring  LI  during  peak  eastward  tidal 
flow.  The  lee-wave  release  mechanism  proposed  by  Maxworthy  (1979)  is  a  plausible  generation 
mechanism  because  it  implies  that  a  hydraulic  control  develops  over  the  ridge  at  LI  that  prevents 
waves  from  propagating  westward  until  the  tide  becomes  subcritical.  This  mechanism  is  unlikely 
because  waves  appear  to  be  "released"  during  peak  easward  tidal  flow  at  the  ridge  throughout  the 
spring-neap  cycle  during  which  time  the  peak  ebb  currents  range  from  strongly  supercritical  to 
subcritical,  and  the  lee-wave  mechanism  is  not  possible  unless  the  tidal  currents  are  supercritical. 
Zhao  and  Alford  (2006)  proposed  an  alternate  mechanism  in  which  waves  are  generated  over  the 
eastern  ridge  during  peak  flood  tidal  currents  (westward),  and  these  waves  steepen  into  trains  of 
nonlinear  waves  as  they  propagate  into  the  basin.  Although  they  proposed  a  different  mechanism, 
Zhao  and  Alford  (2006)  also  found  that  the  generation  phase  was  independent  of  the  strength  of  the 
tidal  currents. 
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The  theoretical  analysis  that  we  have  developed  is  based  on  the  hypothesis  that  the  generation 
mechanism  is  dominated  by  linear  processes,  and  this  is  based  on  the  evidence  that  wave  arrival  time 
at  a  particular  mooring  in  the  basin  is  roughly  the  same  relative  to  the  phase  in  the  particular  tide  that 
generated  that  wave.  That  is,  diumally-generated  A-waves  (Ramp  et  al.  2004)  arrive  at  mooring  S7  at 
the  same  time  every  day,  while  semidiumally-generated  B-waves  (Ramp  et  al.  2004)  arrive  roughly 
one  hour  later  each  day.  This  implies  that  hydraulic  effects  at  the  generation  site  are  weak  and  that 
nonlinear  effects  on  the  propagation  speed  are  also  weak.  These  justify  development  of  a  linear  model, 
which  is  further  simplified  by  assuming  that  the  internal  tides  are  generated  and  propagate 
predominantly  as  first-mode  waves,  which  is  based  on  evidence,  both  from  modeling  and  observations, 
that  more  than  90%  of  the  internal  tidal  energy  in  the  Luzon  Strait  is  in  the  first  mode  (Simmons 
2009). 

The  linear  model  is  derived  from  the  inviscid  equations  of  motion  with  the  Boussinesq  approximation 
for  a  two-layer  fluid  with  upper-layer  depth  d i  and  lower-layer  depth  d2,  with  the  free-surface  defined 
by  rj  and  the  interface  between  the  layers  defined  by  £  and  with  the  bottom  bathymetry  defined  by  h, 
as  shown  in  Figure  2. 


Figure  2:  Schematic  of  the  two-layer  system  used  to  develop  the  linearized  equations. 


The  equations  of  motion  for  the  velocities  in  the  upper  layer,  ui,vi,  and  lower  layer,  U2,  V2,  along  with 
continuity,  are  given  by 
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where /is  the  Coriolis  parameter  and  we  have  assumed  quasi  two-dimensionality  by  ignored  all 
gradients  of  v  and  all  gradients  with  respect  to  y.  After  linearization  and  assuming  a  rigid  lid,  weak 
topography  ( a=h/d,2«l )  and  a  shallow  upper  layer  ( J3=d]/d2«l ),  the  nondimensional  equation 
govering  the  behavior  of  the  interface  when  subjected  to  a  half-cosined-shaped  ridge  is  given  by  the 
inhomogeneous  Klein-Gordon  equation, 

Cm  +  -  y~^Crr  =  ~Af(x)  sin (kbx)  cos(f)  , 


where  R=co/f  is  the  Rossby  number  and  co  is  the  tidal  frequency,  kb  is  the  nondimensional  wavenumber 
of  the  topography.  Here,  x  is  the  nondimensional  distance  and  t  is  the  nondimensional  time,  such  that, 
if  *  represents  a  dimensional  quantitity,  then  x=kx*  and  t=cot*,  and  k=2n/X  is  the  internal  tidal 
wavenumber  and  A  is  its  wavelength,  for  /3=0  and  R=  go  .  In  the  governing  equation  for  z,  f(x)  =  H(x  - 
7i/kb)H(-x  +n/  kb)  and  H(x)  is  the  heaviside  function,  and  the  amplitude  is  given  by 

A  =  , 

where  in  the  present  analysis  Fo=l.  The  solution  for  /is  given  by  a  standing  wave  in  the  region 
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and  freely-propagating  waves  that  radiate  away  to  the  left  and  right  of  the  ridge  (located  at  x=0), 


Cl  =  Coo°s(fcox-|-t)  (x  <  —n/kb) , 
C r  =  —  (0  cos(A-qj  —  t)  (x  >  n/kb) , 


where  the  amplitude  of  the  waves  that  radiate  away  is  given  by 
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and  the  function  G  is  given  by 
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The  nondimensional  wavenumber  of  the  waves  in  the  presence  of  rotation  is  given  by  the  dispersion 
relation  for  Poincare  waves 
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where  c;  =(l+(3)'  .  This  solution  shows  that  the  phase  in  the  tide  in  which  troughs  in  left-  or  westward- 
propagating  linear  waves  cross  the  sill  is  strong  function  of  ko/kb,  since  the  function  G  is  not 
monotonic.  For  most  problems  of  interest,  ko/kb  is  small  since  the  scale  of  the  topography  is  typically 
of  the  order  of  or  smaller  than  the  internal  tidal  wavelength.  However,  different  values  can  lead  to 
wave  generation  at  different  tidal  phases.  Figure  3  depicts  the  analytical  result  for  four  different  values 
of  ko/kb,  and  shows  that  for  ko/kb=0.5, 1.5,  and  3.5  (Figures  3(a),  (b)  and  (d)),  left-propagating 
depressions  coincide  with  the  end  of  the  ebb  (eastward)  tide,  while  for  ko/kb=2.5  (Figure  3(c)),  left- 
propagating  waves  coincide  with  the  end  of  the  flood  (westward)  tide.  Therefore,  because  the  value  of 
ko/kb  affects  the  sign  of  G,  the  ratio  of  the  topographic  scale  to  the  wavelength  of  the  internal  tide 
determines  the  apparent  phase  in  the  tide  at  which  depressions  are  generated.  Short  topography 
relative  to  the  internal  tidal  wavelength  yields  westward-propagating  waves  that  originate  at  the  end  of 
the  ebb  tide,  while  longer  waves  can  alternate  between  the  end  of  the  ebb  and  the  end  of  the  flood  tide. 

These  results  are  contrary  to  the  observations  because  they  indicate  that  waves  are  generated  during 
slack  rather  than  peak  tides.  We  hypothesize  that  the  linear  theory  will  show  that  a  two-ridge  system 
can  lead  to  left-propagating  waves  that  are  generated  during  peak  ebb  tides  over  the  eastern  ridge,  and 
that  nonlinear  effects  are  not  significant  in  altering  this  behavior.  This  is  the  subject  of  ongoing  work. 
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Figure  3:  Contours  of  £=  -(0,  0.25,  0.5,  0. 75,  0.9  5)  (jo  depicting  the  influence  of  the  relative  width  of 
the  sill  on  the  trajectories  of  depresssion  waves  that  emerge  over  the  course  of  three  tidal  cycles.  The 
dashed  line  indicates  the  direction  and  relative  magnitude  of  the  barotropic  currents  while  the  dark 
solid  lines  trace  out  the  peak  depressions  as  they  propagate  to  the  left  and  right  of  the  sill.  The  small 
plots  indicate  the  shape  of  the  bathymetry  for  reference. 


IMPACT/APPLICATIONS 

Understanding  internal  tides  and  waves  is  important  for  developing  parameterizations  of  their  effects 
in  larger-scale  ocean  models. 
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